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Abstract 

The dynamics of a (nonlinear) Berger plate in the absence of rotational inertia are considered 
with inhomogenous boundary conditions. In our analysis, we consider boundary damping in two 
scenarios: (i) free plate boundary conditions, or (ii) hinged-type boundary conditions. In either 
situation, the nonlinearity gives rise to complicating boundary terms. In the case of free boundary 
conditions we show that well-posedness of finite-energy solutions can be obtained via highly non¬ 
linear boundary dissipation. Additionally, we show the existence of a compact global attractor for 
the dynamics in the presence of hinged-type boundary dissipation (assuming a geometric condition 
on the entire boundary [23]). To obtain the existence of the attractor we explicitly construct the 
absorbing set for the dynamics by employing energy methods that: (i) exploit the structure of the 
Berger nonlinearity, and (ii) utilize sharp trace results for the Euler-Bernoulli plate in [25] . 

We provide a parallel commentary (from a mathematical point of view) to the discussion of 
modeling with Berger versus von Karman nonlinearities: to wit, we describe the derivation of each 
nonlinear dynamics and a discussion of the validity of the Berger approximation. We believe this 
discussion to be of broad value across engineering and applied mathematics communities. 

Keywords-. Global attractors, nonlinear plate equation, well-posedness, asymptotic behavior of 
dynamical systems 
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1 Introduction 

We consider the well-posedness and long-time behavior of solutions to a boundary damped Berger 
plate equation taken in the absence of rotational inertia. Its derivation and relation to the scalar von 
Karman nonlinearity [H nulls] are discussed in depth in Section ITTI We consider the isotropic plate 
to be thin, as is usual in large deflection theory [23j . Thus we consider a two dimensional domain 
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C with smooth boundary F. The partial differential equation (PDE) model of interest is then: 


Utt 


A^i 


7-T. 


Vu ■ Vu dfl 


Au = p(x) 


BC{u) 

it(0)=ito; -«t(0)=Mi. 


in (0, T) X fl 
on r 


( 1 . 1 ) 


The term p € represents static pressure on the top surface of the plate. The physical parameter 

T is fixed and positive, while 7 corresponds to in-plane stretching (7 < 0) or compression (7 > 0). 
In this treatment (i) we consider the case 7 > 0, as is by now standard in treatments of the Berger 
plate [13 [IS], and (ii) we normalize T = 1. The term BC{u) in (11.11) represents the plate’s boundary 
conditions. 

Throughout the paper, we will consider nonlinear damping acting through two distinct types of 
boundary conditions; the nature of damping will be critical to our well-posedness and long-time be¬ 
havior results. 

• Free-clamped boundary conditions with dissipation, denoted (FCD): 


{ Am -I- (1 —= 0 on Fi, 
d^{Au) + {1 - p)B 2 U - piu = \ut\^ut on Fi, 

M = dyU = 0 on Fq, 

where F = Fq U Fi is a disjoint union, Fq 0, and the boundary operators Bi and B 2 are given 
by 

BiU = 2h'iV2Uxy — vfuyy — V 2 UXX = —drrU — (div v) dvU , 


B 2 U — dx [(^1 ^ 2 ) ff ^1^2 i^'^yy 

V = {vi, V 2 ) is the unit outer normal to F, t = {— 1 ^ 2 , is the choice of unit tangent vector along 
F. The parameter pi is nonnegative0. The constant 0 < ^ < 1 has the meaning of the Poisson 
modulus. 

• Hinged boundary conditions with dissipation, denoted (HD): 

M = 0; Am = —D{dvUt) on F, (1.2) 

where the damping function D{-) S C^(K) is monotone increasing, and D(0) = 0. 

Remark 1.1. It is clear that, barring issues of elliptic regularity, one can consider various combinations 
of the above boundary conditions. (This is done, for instance, throughout [Il|.) In the analysis of 
(FCD), we consider a disjoint portion of the boundary to be clamped, primarily to simplify aspects of 
our analysis which are not the principle focus, and to avoid technical issues of elliptic regularity. One 
can take the entire boundary to be free, so long as > 0 . 

Remark 1.2. Another configuration (which is of physical interest) involves the so-called hinged-clamped 
boundary conditions. One can consider the damping acting via moments on a portion of the boundary 

^ Since Fq 0, we could take pi = 0 [HI, however, we retain it for the sake of generality. 
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Fi (again, disjoint from Fq— where the boundary is clamped). This configuration, denoted by (HCD), 
is given by: 

{ u = 0; d,jU = 0 on Fn , . 

(F3) 

u = 0; Am = —D{di,Ut) on Fi 

We remark further on the (HCD) configuration in Remark If .81 below. 

For the Berger and von Karman evolutions, taken with homogeneous boundary conditions (of 
clamped or hinged-type), the theory is rather complete. However, in the two inhomogeneous cases 
considered here the results we obtain differ dramatically from those available for other nonlinear plate 
equations (e.g., von Karman, Kirchhoff); this occurs since the advanced theory for these evolution equa¬ 
tions does not apply. In fact, the inhomogeneous boundary conditions are what differentiate present 
analysis of the Berger plate from those available in other available literature (see Section ra . 

We seek to determine the effect of boundary damping on the dynamics, and seek asymptotic behav¬ 
ior results (when possible) under “standard” assumptions on the plate domain and damping functions. 
We also provide a discussion of the validity of the Berger model under these boundary contributions, 
including a review of the pertinent engineering and mathematical literature. By providing a modern 
analysis of the Berger evolution with dissipation acting on the boundary, we attempt to reconcile the 
hurdles in PDF analysis with those described in the engineering literature. Indeed, the engineering 
literature reviewed here provides numerical contrast between the von Karman and Berger models in 
the presence of each of the standard boundary conditions. 

Notation: We will make use of standard Sobolev spaces and denote norms in iF®(fl) {p = 2) 

as II • ||s and || ■ ||o = II • IIl 2 (o). We will use the notation (•, •)n for inner-products in and (•, •)r 

for those in L^{T) (or a subset of F, as indicated by the subscript where necessary). For simplicity, 
norms and inner products written without subscript are taken to be of the appropriate domain 
(e.g., (•, •) on H and (•, •) on F). The space Wq’^(D) denotes the completion of in the kF'*’^(r2) 

norm. We will use the ubiquitous constant C > 0, and denote critical dependencies when necessary. 

1.1 Modeling Considerations 

This portion of the paper is expository. It discusses the origins of the von Karman and Berger non- 
linearities, and provides an engineering perspective on the validity of each model. 

Let u'‘{x,y,z,f) denote the ith component of the displacement vector of the middle surface of 
an elastic plate occupying the region 12 C Let and cr^ denote the standard strain and stress 
tensors, respectivelyj. We assume the plate is homogeneous and isotropic (with associated stress-strain 
relation). 

The so called full von Karman system is based on the employment of a variation principle applied 
to the energy, and the assumption of finite elasticity [231 Ch.2, pp. 13-20]. This is a modification of 
the Kirchhoff approach, which is itself based on (i) a linear strain-displacement relation and (ii) an 
assumption that the linear filaments of the plate remain perpendicular to the central plane of the plate 
throughout deflection. (The second hypothesis is then linearized to produce the Kirchhoff equation.) 
The full von Karman system is arrived at by replacing (i) above with the nonlinear strain-displacement 

^Here, and below, we utilize the summation convention and identify the spatial variables x = xi, y = X 2 , and z = X 3 . 
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relation 


1 / du’‘ du^ \ 1 du^ dvr 

2 \9a;j dxi / ”*" 2 dxi dxj ' 

The resulting system (which has a nice vectorial representation—see |22| and references therein) 
consists of a nonlinear beam equation in the transverse displacement, coupled (nonlinearly) to a sys¬ 
tem of elasticity for the in-plane displacements; much work has been done on this system [22l [29] 
(and references therein). To simplify the structure of these equations, one may neglect the in-plane 
acceleration components of the model [Ml dSj. This allows for the decoupling of the in-plane and 
transverse dynamics. The resulting nonlinearity is referred to as the scalar von Karman nonlinearity: 

fv{u) = -[u,v{u) + Fo], 

where Fq is a given in-plane load. The von Karman bracket is employed: 

[ri, u] = UxxVyy VyyVxx ‘^VxyVxy (^■^) 

for all functions u,v sufficiently smooth [14]. The Airy stress function, v{u) solves 

A'^v = —[u,u] in ri; v = d^v = Q on T. 

Remark 1.3. We note that in the studies of nonlinear plates often rotational inertia in the filaments 
of the plate (represented by the term —{h?'/Vl)Autt on the LHS of the plate equation (II.IL where h 
is the thickness of the plate) is neglected. When h is small (the plate is thin), this term is discardable 
[M] [23]. Mathematically, however, this term is highly non-trivial and has substantial bearing on 
well-posedness and long-time behavior results. 

Even with the simplihcation to a scalar equation, the complex structure of the scalar von Karman 
nonlinearity (making use of a nonlinear elliptic solver) provides incentive to attempt further simplifi¬ 
cation. We now paraphrase a review in [28] : 

In 1955, Berger proposed in m a modified von Karman system. This simplification is based upon the 
assumption that the so called second strain invariant of the middle surface is negligible. The resulting 
equations are simpler than von Karman’s and Berger compared his results to known solutions, and 
found good agreement. This approach was generally accepted, though the hypothesis lacked a clear 
mechanical interpretation. Yet, in many problems solved via the Berger method the edges of the plates 
were assumed to be restrained from in-plane movements, and, numerically, the exactness of the method 
is associated to such restrictions. For freely movable edges the aecuracy of the results furnished by the 
Berger approach becomes questionable. Indeed, the resulting [numerical] accuracy depends on the order 
of the differential operators appearing in the boundary conditions. In this way, hinged and free boundary 
conditions provide more opportunities [than clamped conditions] for the solutions to become inaccurate 
or singular. 

Another analysis in |32j attempts to hnd a rational mechanical basis for the Berger method. How¬ 
ever, the authors also conclude that the Berger results may not be a consistently valid approximations 
across types of boundary conditions. In [5], it is elaborated that Berger’s line of thought leads to 
questionable results for movable edge conditions owing to the fact that the neglect of second strain 
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invariant with movable edges fails to imply freedom of rotation in the middle plane where the elastic 
stress exists. The novelty of [5] is that the authors proceed to provide an additional nonlinear term 
in the equations of motion to correct this, which they report has good accuracy with known solutions 
for all types of boundary conditions. 

1.2 Previous Results and Mathematical Motivation 

We now provide a general overview of previous mathematical work done on the system in (11.11) . 
Note that none of the results below focus on inhomogenous boundary conditions] to our knowledge 
no previous works have provided a PDE analysis of the Berger plate with dissipation acting on the 
boundary—though such analyses are prevalent for the more “complex” von Karman evolution. 

The classical work in [16] establishes results on the existence of exponential attractors for a more 
general Berger-like equation in R” (which subsumes the standard Berger model), taken with homoge¬ 
neous hinged boundary conditions and fully-supported interior damping. This reference discusses early 
work done by Ball (see references in [16j l on extensible beam equations. The methods utilized by the 
authors of m rely on the development of modified energy functionals which allow one to circumvent 
the need for an explicit construction of a Lyapunov functionl^ 

In the more recent book [ID, the Berger evolution with homogeneous hinged boundary conditions 
(along with more general nonlinear terms) is considered with fully-supported interior, linear damping. 
Utilizing a Lyapunov approach, the existence of smooth attractors is shown, but only in the presence 
of (regularizing) rotational inertia terms. 

In |7] the authors study a non-dissipative von Karman model arising in the context of piston theory 
EKIl] for nonlinear flow-plate interactions. This work also considers the case where rotational inertia is 
absent, and the primary difficulties arise via (i) free boundary conditions, as well as (ii) non-dissipative 
(piston theory) terms appearing (11.11) . The analysis does not concern the Berger evolution, yet the 
authors must also combat non-dissipative terms in the setting of free boundary conditions. 

In [26], a Berger beam is considered with free-clamped boundary conditions; in this work stabi¬ 
lization is accomplished via a nonlinear boundary feedback acting in the third order free boundary 
condition. We also point out that to accommodate the non-dissipative terms which arise due to the 
interaction of the geometric nonlinearity and the free boundary condition (as discussed herein), the 
author incorporates the nonlinear boundary contribution into the boundary condition at the free end 
of the beam. Whether this is physical is unclear, however it does present an interesting mathematical 
consideration which could be investigated for free-clamped Berger plates. 

The analysis in [19j deals with the boundary stabilization of a Kirchhoff (polynomial) nonlinear 
plate, taken with rotational inertia, and linear damping via a hinged boundary condition (analogous 
to what we take). Though the analysis is somewhat simplified (due to rotational inertia), this work 
provides a guide for handling the non-dissipative, second-order trace terms in the case of hinged 
boundary conditions below. 

Finally, in perhaps the most modern account, the expansive monograph m addresses the well- 
posedness and long-time behavior of the scalar von Karman equations. However, this text also contains 
preliminary material on linear plates, and an abstract presentation of nonlinear plates. Owing to this, 

®This is a key issue below: in many nonlinear plate/extensible beam problems the construction of a (strict) Lyapunov 
function requires a “unique continuation” result, which may not be generally available [13 HE]. 
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many of the results on well-posedness and long-time behavior are valid for Berger’s evolution in the 
case of clamped or hinged type boundary conditions. However, if any type of free boundary condition 


is employed (homogenous or inhomogeneous), the analysis in [14], although vast, does not apply. 


Moreover, the best results in [14j for the von Karman evolution (with no rotational inertia terms) with 


hinged dissipation (HD) concerns local compact attractors. (In this reference they actually address 
the more general (HCD) conditions—allowing for a clamped portion and hinged dissipation portion 


fRemark ll.2l) —with a geometric condition in force on Fq; their results apply to (HD) with no geometric 


condition by simply taking Fg = 0.) To obtain attractors in a global sense, additional assumptions 
must be made (fully-supported interior damping), and the resulting approach is indirect. For a more 
detailed discussion relating these results to those here, see the discussion following the presentation of 
Theorem 11.31 below. 

In this treatment we are analyzing two cases where the Berger dynamics have a marked differ¬ 
ence from the von Karman dynamics. In the case of free-clamped (FCD) type boundary conditions, 
Berger’s nonlinearity produces an additional non-dissipative term in the dynamics which can only be 
accommodated via a highly nonlinear damping to obtain well-posedness of the model. As such. Section 
[3| can be viewed as a quantifiable mathematical manifestation of the difficulties described above in 
applying Berger’s approximation to the von Karman model in the presence of a free boundary com¬ 
ponent. For our analysis, we assume Berger’s hypothesis is in force and we analyze (ED directly from 
the equations, yielding non-dissipative effects which enter the analysis at the level of finite energy. 

In the case of hinged-type boundary conditions (HD), after showing the well-posedness of the model 
(similar to what is done for (FCD)), we are interested in the asymptotic-in-time analysis. The standard 
nonlinear dissipation acting via moments on the boundary is sufficient to obtain ultimate compactness 
of the dynamics (in the form of a global attractor) without any additional damping or non-standard 
assumptions —though we do require a somewhat classical geometric assumption on the domain (albeit 
for a new reason: the interaction between (HD) and the nonlinearity). In fact, we explicitly construct 
the absorbing ball using techniques from [I2j . and along with modern techniques in the theory of 
dynamical systems (related to the asymptotic smoothness property [141 [2T|, see Sec. 4.2), we obtain 
the existence of a compact global attractor. 

1.3 Functional Setting 

The natural energy for linear plate dynamics is given by 



Here, a(u,v) represents the potential energy and given by the bilinear form 



(1.5) 


where 



( 1 . 6 ) 
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Again, the notation [u,v] corresponds to the von Karman bracket (as introduced earlier in the text), 
and for all functions sufficiently smooth [14]: 


f [u,v]dil, = — [ {{Biu)d^v — {B 2 u)v) dV. 

JQ Jr 


Typically /i < 1/2, but we allow for /i < 1 and note that = 1 corresponds to the limiting case for the 
boundary conditions, and the bilinear form a(-, •) collapses to the standard bilinear form (A-, A-)o in 
the case of (HD) boundary conditions. 

The above energy E{t) dictates the state space 


jr = 


X L^{n) 


for (FCD), 


E{t) = Bit) +-\\yuit)\f 


{H^ n H^){n) X L‘^{n) for {hd). 

Throughout this treatment, we will use the following nonlinear energies associated to (HI): 
S’{t) = S’{u, ut) = E{t) + n(M(t)), 

where H represents the non-dissipative and nonlinear portion of the energy 

n(u) = i(||Vu||4-27||Vu||2-4 [ pu). 

4 Jn 

1.4 Main Results and Discussion 


(1.7) 


( 1 . 8 ) 


In this treatment the main results concern: well-posedness under (FCD) and (HD) boundary con¬ 
ditions, as well as long time behavior of the dynamics corresponding to dnD under (HD) boundary 
conditions. 

Firstly, consider equation dH) with (FCD). We now give a novel well-posedness result (see Section 
[2|for formal definitions of solutions): 

Theorem 1.1. With reference to HI), taken with (FCD) boundary conditions and with initial data 
(uo,ui) € M’, for any T > 0 there exists a unique generalized solution u G C{0,T] J^f) satisfying the 
energy equality 


'S{t)+ J j |Mt|®-l- J J- \\Vu\f){d^u)ut = (F{s) 


(1.9) 


s Ti 


s Ti 


for every t > s > 0. The dynamics depend continuously on the initial data. This implies that the map 
(u(0), Ut(0)) i-G (u(t),ut(t)) defines a strongly continuous semiflow SF{t) on 

Remark 1.4. We emphasize that with only linear boundary damping the global existence of the solutions 
is an open problem, as one cannot obtain the necessary a priori bounds. To overcome this difficulty 
we consider a nonlinear boundary dissipation which has an extreme (yet minimal with respect to our 
analysis) polynomial structure. 

Secondly, we analyze the Berger model with (HD) boundary conditions. It should be noted that 
under (HD) conditions (unlike the (FCD) conditions) the structure of equation (11.11) does not yield any 
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non-dissipative terms in the energy relation. This makes the direct application of the abstract plate 
theory more straightforward for obtaining a well-posedness result. Before giving this, the following 
assumption is needed on the damping: 

Assumption 1. There exists positive constants 0 < m < M < oo such that 

m < D'{s) < M, |s| > 1 

Remark 1.5. For the purposes of well-posedness alone, one could impose weaker assumptions on the 
damping mechanism D{s)] e.g., D{s) could exhibit arbitrary polynomial growth. However, for our sub¬ 
sequent results concerning long time behavior it is necessary that D{s) satisfy the stronger assumption 
above which is also utilized in nnn]. 

Now we give the well-posedness result of dm taken with (HD) boundary conditions: 

Theorem 1.2. Let As sump tion\^ hold. With reference to il.ll) taken with (HD) boundary conditions 
and with initial data (mo,Mi) G , for all T > 0 there exists a unique generalized solution u G 
C{0,T]J^) depending continuously on the initial data. This implies that the map (m( 0), itt(O)) i—>■ 
{u{t),ut{t)) defines a strongly continuous semiflow Snit) on J^. The following energy equality holds: 

t 

J J D{di,ut) {dj,ut) = S' (s). ( 1 - 10 ) 

s r 

Remark 1.6. Assumption [T] provides linear bounds for the damping function from below and above 
which ensures the validity of the energy equality (for more detail, see [Ill Section 4.2]). 

Existence and properties of compact attractors for nonlinear plate models (e.g., von Karman, 
Berger, Kirchhoff) have been studied in various contexts over the last 30 years0 In the present case, 
the dynamical system associated to the Berger dynamics with (FCD) cannot easily be shown to be 
dissipative. The absence of a strict Lyapunov function (not readily available) seems to preclude showing 
the existence of a compact global attractor for the dynamics. However, we can show dissipativity of 
the dynamical system under (HD) boundary conditions. Thus, by showing the asymptotic smoothness 
property of the associated dynamical system, the plate dynamics under (HD) will admit the existence 
of a compact global attractor for finite energy, generalized solutions. For this, we need an additional, 
standard, star-shaped geometric assumption on the boundary mug: 

Assumption 2. There exists a point Xq G such that h(x) = x — Xq has the property that h ■ v > 0 
for all X G F, where v is the outward normal vector to F. 

Remark 1.7. The star-shaped condition given above arises in boundary control, when the boundary 
(or just a portion) is subjected to a given feedback. In using state-of-the-art trace estimates [25] . 
we can accommodate boundary terms associated with the linear dynamics without a geometric condi¬ 
tion. However, to deal with the nonlinear Berger contribution on the boundary (in utilizing the flux 
multiplier h ■ Vit), we utilize the star-shaped assumption on F. Throughout [23] a star-shaped and 
star-complemented geometric condition (what is in Assumption |2| with the additional assumption that 

^See Section o also see the Appendix for a discussion of dynamical systems and concepts relevant to the treatment 
here. 
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h['x.)-v < 0 on To) is in force for obtaining uniform decay rates. Such a star-shaped, star-complemented 
assumption is also used for various configurations in the context of von Karman equations in |14) . We 
note that our assumption on T is precisely what [25) can eliminate for linear dynamics while still 
obtaining the same results of [23] : however, at present it seems indispensable to have an assumption 
on the boundary F due to the interaction of the nonlinearity and the boundary conditions. 

We note that at present the best result obtained for non-rotational von Karman dynamics in the 
(HD) configuration (with additional damping and assumptions) is that of local attractors 1 |141 Theorem 
10.5.7]). This is to say that when one restricts the semiflow to a ball of radius R of initial 

data, a compact attractor can be shown for the restricted dynamical system. This attractor is not 
necessarily uniform with respect to large R. Such a result follows (though not immediately) from the 
so called asymptotic smoothness property of the dynamics, and estimates which yield a local absorbing 
ball. (See [T31 Section 10.5] for this presentation.) To obtain a global attractor the authors of [H] utilize 
the fact that a gradient dynamical system, which is also asymptotically smooth, has a compact global 
attractor fTheorem 16.21 in the Appendix here). To show the gradient property of the dynamics, they 
rely upon (additional) fully-supported interior damping. This damping yields a “unique continuation” 
property for the dynamics (or, from another point of view, allows for “uniqueness” in a compactness- 
uniqueness argument). This approach, and the requisite additional damping, is necessary because 
“direct multipliers” do not yield a uniform absorbing ball for the von Karman dynamics—owing to 
boundary contributions and the structure of the H for von Karman dynamics. 

In the case of Berger’s nonlinearity, we are able to adapt an approach in [12] to show the existence 
of an absorbing set for the entire dynamics. This is precisely because of the structure of the Berger 
nonlinearity. The discussion above brings us to the primary result concerning long-time behavior of 
the Berger dynamics with (HD) boundary conditions: 

Theorem 1.3. Consider taken with (HD) boundary conditions; under AssumptionsUl and\^ the 
dynamical system (Snit), J^) has a compact global attractor. 

Remark 1.8. The above holds in the case of (HCD) boundary conditions (see Remark B.2L where 
r = To U Fi and the dissipation is active only on Fi. This will be addressed in the forthcoming 
manuscript [I] , which will also investigate further properties (such as dimensionality and smoothness) 
of the attractor given in Theorem 11.31 The key point in analyzing the partially-damped boundary 
configuration is to address the higher-order trace term Au on the uncontrolled portion of the boundary 
Fq. a geometric condition on the entire boundary (for instance, the star-complemented, star-shaped 
condition mentioned in Remark ini) is not sufficient in obtaining estimates which control the size of 
the absorbing ball. Rather, precise control of Aujro must be obtained via techniques which can be 
found, for instance, in |2|3] (and references therein). 

1.5 Technical Challenges and Contribntions 

With regard to much of the multiplier analysis, Berger’s nonlinearity behaves more nicely than von 
Karman. And, in the case of homogeneous hinged or clamped boundary conditions, the analysis of the 
von Karman plate essentially subsumes the analysis of the Berger evolution. However, under (FCD) 
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boundary conditions (with Fq G Hq^^I)), we have the disparity: 


{fv{u),ut)n = (-[m,u(m) + Fo],Mt)n = ~ |i||Au(M)||^ - ([u,Fo],m)o| 
{fB{u),ut)n = ((7- ||Vu||^)Au,Mt)j^ = ~ |i||Vu||^ -7||Vu||^| 


(7- ||Vu|n(aj.M,ut)ri 


( 1 . 11 ) 


( 1 . 12 ) 


The boxed term highlights the fundamental difference of the dynamics, as the Berger-free dynamics 
contain an inherently non-dissipative term with which we must contend. 

In the case of (HD) boundary conditions, we note that the terms of interest arise from the bi- 
harmoiric term—which plays a key role in the multiplier analysis. We have for any tp sufficiently 
smooth 

(A^m, (j))Q = (Am, A(j))n + {d^{Au), (j))r + {D{d^ut), d^(j))r (1-13) 

In this case oi cp = h-'Vu (for h a “nice” vector field, u G H^(fl)) the damping on the boundary is pitted 
against a high order boundary term, which must be accommodated for any sort of long-time behavior 
analysis. Such trace terms must be handled delicately. Here, we are able to (via the special structure 
of the Berger nonlinearity) control these boundary terms directly, in a way that allows techniques from 
[T^ to obtain. 


1.5.1 Principal Challenges 

The non-dissipative term above for (FCD) in (jl.l2|) appears in the energy relation for the dynamics 
(and it is absent for von Karman dynamics). We follow the modern approach in |14) . using the abstract 
plate theory developed therein, and we must utilize fully nonlinear damping of high polynomial degree. 
As such, we must address the technical difficulties associated to such damping when utilizing energy 
methods. For obtaining dissipativity and asymptotic smoothness of the dynamical system associated 
to the (HD) dynamics, we need to handle the nonlinear damping carefully. 

In both configurations—(FCD) and (HD)—obtaining the results in Section [F4l requires addressing 
the criticality (with respect to the finite energy topology) of the Berger nonlinearity. Similar to the 
non-rotational von Karman dynamics, the Berger nonlinearity is not compact, and this requires some 
care in energy estimates. 

Finally, what is perhaps the most challenging aspect of the analysis herein, is the difficulty involved 
in handling the trace terms associated to (HD) dynamics in the long-time behavior analysis. These 
terms include the standard linear terms—due to the principal biharmonic part—as well as “new” 
boundary terms arising due to the interaction of the boundary conditions, requisite multipliers, and 
the Berger nonlinearity. We must deal carefully with signed terms on the boundary (including imposing 
a geometric assumption on the domain), as well as utilize the sharp trace results on second order traces 
(as in [25]). 

1.5.2 Principal Contributions 

In what follows we address each of the points outlined above, as well as make additional contributions: 
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• We show that a high (polynomial) degree of nonlinear damping can accommodate (at least in 
terms of well-posedness) the free-type boundary condition (FCD). 

• We demonstrate asymptotic smoothness for the dynamical system associated with the non- 
rotational Berger dynamics under (HD) conditions. This utilizes a decomposition of Berger’s 
nonlinearity, akin to the one utilized in the analogous result for non-rotational von Karman 
dynamics [TU pp. 496-497]. 

• We critically use what is in m in considering (HD) conditions, but we have a different treatment 
of trace terms (particularly in the use of the sharp estimates in [25]). This has the benefit of 
extending the approach in m to the non-rotational plates considered here. 

• For the (HD) case, we are able to adapt the approach in [12] (corresponding to a semilinear 
wave equation with boundary damping) to construct the compact global attractor by directly 
obtaining estimates on the absorbing ball. This depends in a critical way on the special structure 
of the Berger nonlinearity (and does not seem possible for von Karman). We point out that, in 
some sense, we obtain a stronger result than in the analogous case for von Karman in M- 

Remark 1.9. The authors of [14] utilize additional assumptions on the damping to obtain a 
unique continuation result that provides a strict Lyapunov function for the dynamics. Then, 
they use an indirect result (for a gradient dynamical system, asymptotic smoothness implies the 
existence of a compact global attractor—see the Appendix, Theorem 16.1|1 : this is to circumvent 
the need for direct estimation on the absorbing ball in the case of non-rotational von Karman 
dynamics. 

• Most importantly, by addressing non-rotational Berger dynamics with boundary dissipation we 
provide results which seem to be conspicuously missing in the mathematical literature. As we 
have pointed out, the case of (HD) and (FCD) seem to be pathological, in that the analysis of 
the Berger dynamics is not subsumed by that of the von Karman. 

2 Abstract Plate Equations and Well-Posedness 

The issue of well-posedness of strong or hnite energy solutions for nonlinear plates is recent. For a 

detailed and complete discussion of general, abstract plate equations the reader is referred to [T^ . 

For clarity, let us recall some definitions and results. The general model arising in nonlinear plate 


dynamics is: 



( 2 . 1 ) 


To address the well-posedness for problem m, we must give assumptions and properties of the 
operators. 

Assumption 3. With reference to problem 12.1\) : 

• Let A be a closed, linear, positive, self-adjoint operator acting on a Hilbert space %, with ^{A) C 


n. 
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• Let U be another Hilbert space, and Uq be a reflexive Banach space, such that Uq C U Q Uq. 
Additionally, assume that y : Uq ^ Uq is a continuous mapping such that g(0) = 0 and 

(g{vi) - g{v 2 ), Vi - V 2 ) > 0, for all vi,V 2 S Uo 

where {■, ■) denotes the scalar product on U, or the duality pairing between Uo and Uq. 

• The linear operator G : Uq ^ % has that \ Uq % is bounded (or equivalently, 

G*A : Uq is bounded, where the adjoint operator G* is defined by the relation 

{G*u, v) = (m, Gv) ). 

• Let V be a Hilbert space satisfying GV dTL gV G with all injections being 

continuous and dense. The nonlinear operator F : x V ^ V is locally Lipschitz; that 

is: 

\\F{ui,vi) - F{u 2 ,V 2 )\\v' < L{K) ( A^^'^{ui-U 2 ) + ||wi --^211 ^) 

^ 'hL ' 

for all {ui,Vi) G 3l[A}-l'^) x V such that , ||?;i||y < K. 

Remark 2.1. In concrete applications, the term AG'g{G*Aut) models the boundary dissipation. Here, 
the operator G typically represents a suitable Green’s map, and the mapping g : Uq ^ Uq is a. 
Nemytskij-type operator (see Section 3 and 4). 

If we rewrite problem (12.11) as a first order equation we have the equivalent problem 


—u{t) + Au{t) = 


0 


u(0) = Uo = {uq,ui), 


L F{u(t),ut{t)) 

where u(t) = {u{t),ut{t)) and the operator A : ^{A^/'^) xV^ x H is defined by 


r 0 -I 
~ I A AGgG*A \ ’ 

where ^(A) consists of elements (x, y) G x possessing the property 

A{x V Gg{G*Ay)) G v'. 


( 2 . 2 ) 


This structure of problem (12.11) leads the following definitions of strong and generalized solutions in 
this setting. 

Definition 1. A function u(t) is said to be a strong solution to 12.1\} on a semi-interval [0,r), iff 

. u(t) GC([o,r);^(^i/2))ncH[o,r);H) 

• u G b; jA^A}^"^)) and Ut G b; V) for any 0 < a < b < T 

• A[u{t) + Gp{G*Aut{t))] G V for almost all t G [0,r] 

• Equation (12.11) is satisfied in V for almost all t G [0, T] 

• The initial condition in mi holds. 
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The function u(t) is a strong solution on [0,T] if in addition, we have that {u{t),ut(t)) is continuous 
at t = T. A function u{t) is said to be a generalized solution to (12.111 on [0,T], iff 

• u{t) € C(0, T; ^(.Ai/2)) n C\0, T; V) 

• The initial condition in HU holds. 

• There exists sequences of strong solutions {un{t)} to problem \2.1]) defined on [0,T] with initial 
data (uon,uin) such that 


lim max 1 \\unt{t) 
n—^oo tg[o,T] I 


Uti 


ly ■ 


A^^'^{Un{t) - u{t)) I 

n) 


= 0 


The function u{t) is the generalized solution on a semi-interval [0,T), if u(t) is a generalized 
solution on each subinterval [0,r] C [0,r). 


The following theorem provides the existence and uniqueness of local and global solutions. This is 
the primary result of Section 2.4 in [14]. 

Theorem 2.1. Under Assumptions^ the following hold: 


• Local strong solutions: For every uo,ui € such that A[uq + Gg{G*Aui)] G V there 

exists fniax > 0 and unique strong solution such that {u,ut) € C'([0,tmax) X V). 

• Local generalized solutions: If {uo,ui) G ^(A) C x V, where ^(A) is the closure 

of &{A) in X V. Then there exists tmax > 0 and unique generalized solution such that 

eC([0,t„,ax);^(.4i/2) X V). 

• Global solutions: If in addition, strong (or generalized) solutions satisfy 


sup 




n 




for every existence semi-interval then the local solutions referred to above are global, which 

%S to say ^niax ^ oo. 


3 Berger Model with (FCD) Conditions 

3.1 Well-posedness—Proof of Theorem II.Il 

We consider problem (HU with (FCD) boundary conditions. The proof of Theorem 11.11 is based on 
the application of the abstract setup above and Theorem 12.11 


3.1.1 Application of Abstract Setup and Local Well-posedness 

Firstly, we need to modify the abstract model given above to fit our case. For the functional setup, 
we introduce the following spaces and operators: 

• n =V=L^{n), U = L\Ti) X 
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• Uq = X i^3/2(^l), t/' = ij-i/2(ri) X i/-3/2(ri), 

• Au = A^u, u G ^(^), where 

{ u = Vit = 0 on Fq 
M G + (1 — h)B 2 U — fiiu = 0 on Fi > 

Am + (1 — p)Biu = 0 on Fi 


• V = L^{il) and which is given in Section 1.3. 

• F{u,v) =-fsiu)+p, with /b(u) = (7 - ||Vu||^)Am, 

• G : Uq ^ y. , G{gi,g2) = Gi(gi) — G2{g2) where Gi, i = 1 , 2 , are defined by 


Gi(m) = m iff < 


A^m = 0 in 11, M = Vm = 0 on Fq 
Au + (1 — p)Biu = M on Fi 
di^Au + (1 — g)B2U = piu on Fi, 


and 


{ A^u = 0 in 11, M = Vu = 0 on Fq 
Au + (1 — fi)Biu = 0 on Fi 
d^Au + (1 — g,)B 2 U = piu + m on Fi. 


• The mapping g : Uq ^ Uq has the form 


g ■■ {vi,V2) {giivi),g2{v2)), {vi,V2) G Uo 


With the above notation, taking into account that in our case: gi = 0 and g 2 {v) = G{gi,g 2 ) = 

—G 2 {g 2 )- We reduce the equation in (11.11) to the following form 

utt + .4 [u - G 2 ([|ut|®ut] IrJ] = F{u{t),ut{t)), t > 0, 

or equivalently, 

Utt + A[u + Ctg{G*Aut)] = F{u{t),ut{t)), t > 0. (3.1) 

Now, before giving the well-posedness result for (EU, we verify that the conditions in Assumption |3] 
are satisfied. 

Lemma 3.1. With reference to 15*. Ill , operators A, G and F introdueed above eomply with the require¬ 
ments of Assumptions^ 

Proof. By definition of A as above, A is closed, linear, positive, self-adjoint, densely defined on TL 
with S!{A}I'^) = iFp^(ll), and we have dense, continuous injections: 

<zV = 'H = V' c 9 {A^/'^)' = [H^^in)]'. 


The Nemytskij operator 

( 51 , 32 ) = ( 0 , 32 ) ■■Uo^U'o 
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is monotone [HllST] and continuous, since C C(ri) by the Sobolev embeddings. 

The requirement on G in Assumption |3] reduces the fact that 


:n^Uo = X (3.2) 

is bounded. An application of Green’s theorem with the calculations given in [Ml Section 3.2.2] yield 
that the mappings 

Gi : G2 : i^(ri) ^ 

GtM : ^ G*2A : 

are bounded and for every u S C We then have that 

G*Au = {GIAu,-G*2Au) = (5^M|ri,n|rJ. 

By the last equality, (E21) translates into the fact that 

u —>■ G*A^^^u = {d^A~^^^u,A~^^^u) 

is bounded from into x This follows from the identification = 

and the trace theorem, which imply that the maps 

5,(-) : H^{n) ^ Hi/2(r), (OIr : ^ 


are bounded. The validity of requirement on G is thus observed. 

Now, we need to verify that the operator F : x F —>• V'{= V), defined by 

F{u,v) = -/b(w) +p, 

is locally Lipschitz. Let u^uu G with ||A^/^u||-h, ||A^/^w||^ < R, where the bilinear form 

a(-, •) plays a special role for the norm in under (FCD) and is equivalent to ||A • ||o for (HD)E1: 


||/s(w) - /bMIIv < G||m - w||®(^i/ 2 ) + 11 llVwlpAu - llVwlpAwl |. 


Now, 


11 ||Vu|pAu — IjVralpAwl I <| | jlVwlpAu — ||Vu|pAty + ||Vit|pAw — ||Vr(;|pAw| | 

<11 ||Vm||^A(u- w) + AwdjVitjl^ - ||Vm;||^)|| 

< G|H|?,a||w - Hb.o + G|H| 2 ,o(||Vu|| + ||Vn;||) (||Vu|| - ||Vn;||) 

< C'l|u||?,al|w - t«||2,n + G||w||2,n(||Vu|| + ||Vw||) (||m - u;||i,n) 

< G(i?)||u-u;||^(^i/2), 

where we have freely used the Poincare’s inequality, the equivalence F[^{n) = and Sobolev 

^For (HD), this follows from the well-posedness of the biharmonic problem associated to the these boundary conditions. 
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embeddings. The locally Lipschitz property for F is thus observed. 

Moreover, for initial data (uo,ui) £ since the set 

£ = n X 

is dense in and belongs to the domain of corresponding operator A given in (12.211 . the application 
of Theorem 12.II gives the local (in time) existence and uniqueness of generalized solutions. □ 

A Priori Bound and Global Well-posedness—Completion of proof of Theorem \l.l\ 

Having obtained local generalized solutions, our aim is to obtain global a priori estimates for these 
solutions. We note that our computations are first performed on strong solutions, and then extended 
via density to Jff in the limit. Our starting point is the energy relation (lElD, which follows via 
integration by parts for strong solutions and remains valid for generalized solutions. It can easily be 
shown m that there exists positive constants Cq , Ci and C such that 


coS’{u, Ut) — C < E{u, Ut) < ci(a{u, Ut) + C (3.3) 

for all generalized solutions {u,ut) £ JifE Now, using (13.31) in the energy identity (ll.9|) . we have 

t t 

E{t) + j J |ut|® + J J{-f-\\Vuf){d,u)ut<EiO) + C 


0 Ti 


0 Ti 


Dealing with the non-dissipative term above requires some calculation. By the Holder inequality 
we have 


Z 

J y(7-l|Vuf)(a^u)(Mt) 


0 Ti 


I 0 0 

Now, applying the trace moment inequality for d,jU^ 0 Theorem 1.6.6] we have 
t r t t 

J \\^uf){duu)ut <c\j-f WuWy^ llVuiy^ + J llVuf hth^T^) \ ■ 


0 r. 


We note the equivalence of a(-, •) and || • || 2 - We then use Young’s inequality for the first integral on 
the RHS, with powers 2, 4, 8, and 8, so (1/2) + (1/4) + (1/8) + (1/8) = 1. For the second integral of 
RHS, with powers 4, 8/5, 8, where (1/4) + (5/8) + (1/8) = 1, we obtain 


t t I 

J Jil-\\'^uf)id,.u)ut <Cej‘^t + Ce J E{s)ds + e J ||ut||®2(r^) . 

0 Ti 0 0 


This yields 

®For more details, see the discussion preceding ga. 
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t 

E{t) <C{-r,t,E{0)) + C j E{s)ds. 

0 


An application of Gronwall’s inequality gives: 


E{t)<C{T), for all t G [0,r]. 

This implies the global existence of strong (and thus generalized) solutions to equation (j3.1ll . Since 
these solntions also solve the original problem— ( 1131 ) with (FCD) boundary conditions—the proof of 
Theorem 1 1.1 1 is now completed. 

Remark 3.1. At the present time, we do not see how the global well-posedness of solutions can be 
established in the absence of such highly nonlinear damping. As discussed above, this is an indication 
that the dynamical systems properties of solutions does not follow readily from the equations via 
standard techniques used at the finite-energy level. Without highly nonlinear damping active on the 
boundary, we may not have control of the non-dissipative term 

f (7-||Vn||2) [ {d,u)ut drdt. 

Jo JFi 

Additionally, since we have no control over E(t) as T —> -|-oo, long-time behavior analysis is not 
viable. 


4 Berger Model with (HD) Conditions 

In this section we consider problem dnj with (HD) boundary conditions. 

4.1 Well-posedness—Proof of Theorem 11.21 

Similar to Section [3l the proof of Theorem 11.21 is based on the abstract setup given in Section [2] and 
the application of Theorem 12.11 We adapt the abstract model (EU to this case. 

4.1.1 Application of Abstract Setup and Verification of Hypotheses 

We introduce the following spaces and operators: 

• n=L^n), U = L^T) 

• t/o = iJi/2(r) u'o = 

• Au = A^u, u G ^(A), where 

S!{A) = in (iJ^ n ilg)(D) : Art = 0 on T | 

• V' = V = L^(Q) and ^(A^A) = (ffi n E^)(n). 
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F{u, v) = —/b(m) + p where fsiu) = (7 — || Vm||^)Au. 


• G : L^(r) —i> L^{n) denotes a biharmonic extension of the boundary values defined on F. That 
is: 


A^u = 0 in n, 
It = 0 on r, 
Am = m on F. 


G{v) = u iff < 

The mapping g : t/g —>■ Gg is determined by the function £>(•) according to the formula 


= / D{v)wdr v,wGUo. 
’ ° Jr 


We note that the application of Green’s formula gives that G*A = ^|p, and hence ker[G*yl] fl 
contains Hq{Q) and is dense in L^(G). This implies that for all initial data (mo,mi) € 
the set ^{A) x (ker[G*^] n is dense in ,Af and belongs to the domain of the corresponding 

operator A. 

In line with the above setup, the verification of the requirements of Assumption follows as the 
case of (FCD) damping. For more details the reader is referred to [HI Section 4.2]. The proof of 
Theorem 1 1.2 1 is obtained by the application of Theorem 12. II via the density argument above. 


4.2 Long Time Behavior Under (HD)—Proof of Theorem 11.31 

This section is devoted to study the long time behavior of dynamical system generated by the solution of 
cu taken with (HD) boundary conditions. Our principal goal is to prove Theorem If.31 the existence 
of a compact global attractor for the dynamics. 

We utilize the notations consistent with the case of (FCD), but now specified to the case of (HD): 
E{t) = -{llwtlP + ||Au|p}, as well as the identical definitions for E and S' from before. 

As in the proof of [TH Lemma 1.5.4], for all {u,ut) G JS" = (il^ n Hq){V,) x T^(D), and any e > 0, 
it is straightforward to show the potential energy bound (as we will refer to it below): 


|||VM(t)|p + (p,m) 


< e 


llAM(t)lp + illVM(t)l| 


+ M{e,j,p). 


Remark 4.1. In fact, as in m, the key bound which implies (14.11) is 

1 


m\2-ri < e 


||AM(t)ll^ + -l|VM(t)|r 


+ M(e), ? 7 , e > 0 . 


(4.1) 


(4.2) 


This also yields the energetic bounds 


coA(m, ut) — G < S{u, ut) < ciE(u, ut) + G, (4.3) 

for some cg, ci, G > 0 depending on p and 7. Accordingly, we introduce more notation for the study of 
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long-time behavior (following [12)1: 


0 < <gM{t) = S{t) -b M, (4.4) 

where M = is the constant given in (14.1|) . 

4.2.1 Overview of the Proof 

Proof of Theorem 11.31 (which gives the existence of a compact global attractor for {Sh,J^)) follows 
from the application of the abstract Theorem 16.11 fin the Appendix). Utilizing Theorem 16.11 requires 
showing that the dynamical system generated by solutions to (HU) taken with (HD) conditions is 
dissipative (i.e., possess a bounded absorbing set), and has the asymptotic smoothness property. Thus, 
to prove Theorem If.31 we will first present the dissipativity analysis of (11.11) in Section [4.2.21 We will 
utilize energy methods to show the existence of the absorbing ball that attracts every trajectory after 
a certain time. This will be given as Theorem 14.II Then, in Theorem 14.41 the asymptotic smoothness 
property of the dynamical system is given. Its proof follows through several steps in Section 14.2.31 
With these two results, the existence of a compact global attractor for {Sh, will be demonstrated. 

Remark 4.2. We will show an explicit bound on the size of the absorbing ball. We note that this is 
in contrast to the approach taken in m for the non-rotational von Karman plate with hinged-type 
boundary dissipation; there, the dynamics (under additional assumptions, including requisite frictional 
damping in the interior) are shown to be gradient, and thus the existence of a compact global attractor 
(roughly) follows from the asymptotic smoothness property (see Theorem 16.21 in the Appendix). Also 
see the previous discussion following Theorem 11.31 


4.2.2 Explicit Estimates on Absorbing Ball 

In this section, as the first step of the proof of Theorem ll.31 we show the dissipativity of the dynamical 
system corresponding to generalized solutions of (HU with (HD) boundary conditions. To this end, 
we give an explicit estimate on the absorbing ball for which we utilize multiplier techniques and sharp 
trace results for the linear plate equation with a given RHS [25] . We will show the following theorem: 

Theorem 4.1. Let Assumption^ and\^ hold. Then there exists an absorbing set C for general¬ 
ized solutions to HU taken with (HD) boundary conditions. This is to say: for all Rq > 0 and initial 
data (uo,ui) G with ||(ito, ui)||.^ < Rq, there exists a to = t{Ro) such that {u{t),ut{t)) G for 
t>to. 

In order to prove this central theorem, the main technical ingredient will be the following “observabil¬ 
ity” estimate, stated below. 

Lemma 4.2. Let Assumption\^ is in force. Let T > 0 and 0 < a < T/2. Then any solution to (11.11) 
taken with (HD) boundary conditions satisfies the following estimate: 


(T - 2a)E{T) < CE{Q) + C{T)\\D{d,ut)\\l.f^o,T-L^^T)) / {E{T)fdT 


'-JO 


+ C(T) \\D{dyUt)\\i^2^o^T-,L'^(T )) E \\duUt\\i^2^o,T\L'^(T)) EC{T). 


(4.5) 
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Proof of Lemma \4.2[ Now, let m be a generalized solution to (HID taken with (HD) boundary condi¬ 
tions. We begin with the equipartition multiplier u and the flux multiplier h-\7u, where hfx) = x — xq 
for some appropriately chosen xq £ (satisfying Assumption [2|) . Multiplying (11.11) by u and h ■ Vit, 
respectively, integrating over the space-time cylinder, and using Green’s theorem we have 

/ {\\^u\f+ \\'^u\\'^-\\ut\\'^}+[ [ D{d^ut){d^u) =-{ut,u)\^+f {(p, u)-b 7 l|Vu|p} , (4.6) 

Jo Jo Jr Jo 

and 


{||Ait|p-b ||ut|p}-b J y |^|AmP(/i • I^)-b 9^(Au)(h • Vit) - (Am)i9i,(/i • Vit) 

+ 1 f ( 7 -||VM|n [ {h-v)\d^u\'^ =-{ut,h-Vu)\^ + [ {p,h-Vu). 

^ Jo Jr Jo 


(4.7) 


We have used the facts: 


((7- ||Vn|nAM,h-Vu) = (||Vn|p- 7 )(Vm,V(/i-V n)) +y [ 7 -||Vu|p] J {d,u){h-Vu) 

^Vm, V(/i ■ Vu)^ = -div(/i|VMp) and h-Vu = {h ■ v)dvU on F, 

owing to the divergence theorem, and the zero Dirichlet condition on the whole of F (for (HD)). Now, 
taking a suitable combination of (14.6|) and (14.7|) we obtain: 

/ {3||Am||2 -b ||Vu||^ -b llutll^-7l|Vu|p} -b f (7-||Vu|n f {h ■ iy)\d^u\^ 

Jo Jo Jr 

pT pT pT p 

= -{{ut,u)\^+ 2{ut,h-Vu)\^} + {p,u) + 2 {p,h-Vu)- / D{d^ut){d„u) 

Jo Jo Jo Jr 

+ f [ {2(AM)5,y(h • Vm) — (h • r/)|AMp — 2 i9i,(Am)(/i • Vu)}. 

Jo Jr 


(4.8) 


Using the bound on the nonlinear potential energy and reorganizing the terms above, we have a 
preliminary estimate: 


c / E{t) < 
Jo 


(Tl) - ^{ut,u)\^+ 2{ut,h-Vu)\^'^ + J {p,u) + 2j {p,h ■ Vu) + C{j) ■ T 

(T2) -2 [ [ d,{Au){h ■ Vu) - [ {-f -\\Vu\\^) [ {h ■ i^)\d,u\^ 

Jo Jr Jo Jr 

(T3) -[ [ D{d,ut){d,u)+ [ [ 2{Au)d,{h-Vu)- [ [ {h ■ iy)\Au\^ 

Jo Jr Jo Jr Jo Jr 


(4.9) 
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Firstly, by m and Young’s inequality, we observe 


|(T1)| < CE{0) + e 



+ C{p,e,'-f) -T. 


(4.10) 


Now, employing the geometric Assumption [5] on (T2)-(T3), dropping the negatively signed terms 
on the RHS, combining our computations thus far and reorganizing, we obtain the following interme¬ 
diate relation: 


c / E{t) < CE{0) + , 


E{T) + / E{t) 


+ C(p, e, 7 ) • T 


-2 


f f d^{Au){h ■ \/u) - [ ( 7 -||Vm|P) f 
Jo Jr Jo Jr 

[ [ D{d„ut){dj,u)[ f 2{Au)dj,{h -Vu). 
lo Jr Jo Jr 


(4.11) 


In order to estimate the RHS of the above inequality—which involves higher order trace terms—we 
follow the analysis of [19] (which itself critically relies on Theorem l4.3l from [25] given below). However, 
we approach the higher order trace term involving 3y(Au) in a more straightforward way which, 
ultimately, benefits our analysis and fundamentally exploits the structure of the Berger nonlinearity 
to obtain cancellations. 

Now, let the operator A = A, acting on L^(H) with domain i7^(H) fl and let Dl be the 

associated Dirichlet “lift” map defined by 


DLg = f&L^m ^ A/ 


0 in H and f = g on F. 


Let F = I — /b(m) -l-p|. Accounting for boundary conditions, the operator representation of (11.11) is 
then 

utt+A[Au-DL{Au\^)]=F. (4.12) 

Applying A~^ to (|4.12l) (justified on strong solutions, and a posteriori on generalized solutions via the 
corresponding estimate) we obtain 


Au = Au =—A ^utt + Dl (^Au\^^ + A ^F in 

Now, taking the normal derivative of both sides of above equality, multiplying by h ■ Vm, integrating 
over [0, T] x F, and reading off from the equation we have the relation 

f ( dv{Au)[h-Vu) = - f f d^{A~^utt){h-Vu) + ( f O^Dl {Au\^) {h ■ Vu) 

Jo Jr Jo Jr Jo Jr 

+ [ ({d,A-^F){h-Vu). (4.13) 

Jo Jr 

In order to estimate each term of the RHS of ()4.13|1 we use the approach in [19] . The following is the 
critical step which allows our approach to the higher order trace terms to obtain. This relies eritically 
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on the structure of the Berger nonlinearity. If we note that 


A-^F = - -i]^u + A-^p 

= [l|Vu||^ -"i\u + A~'^p, 


(4.14) 


and again h ■ Vm = {h ■ v)dvU^ since it = 0 on F; we have: 


-2 


/o Jt 


d^{Au){h -Vu) - [ ( 7 -||Vu|p) f {h ■ n)\d,^u\'^ 

Jo Jr 

2 / f d,j{A~^utt){h ■ Vu) - 2 f [ (Am|j,) (h ■ Vm) 

Jo Jr Jo Jr 

+ ‘2 f ( 7 -||Vm|P) f {h ■ iy){d^u)'^ - 2 ( f d,y{A~^p{x.)){h-Vu) 
Jo Jr Jo Jr 

- r(7-iiv«i|2) [ ih-,.)\d.u\^ 

Jo Jr 


(4.15) 


Remark 4.3. The cancellation in the nonlinear terms above is critical to the argument below. We will 
be able to use the sign of the term and the standard star-shaped geometric condition Assumption [5] 
on r to discard the nonlinear boundary contribution from our energy approach. 


-2 


to Jr 


= 2 


d,j{Au){h -Vu) - f ( 7 -||Vm|P) f {h ■ i^)\d^u\'^ 

Jo Jr 

( [ d^{A~^utt){h-Vu)-2 f f O^Dl {Au\^) {h-Vu) 

Jo Jr Jo Jr 



-2 / d,y{A ^p{x)){h ■ Vu). 


(4.16) 


Jo Jr 

Now, we estimate the integrals on the RHS of the above inequality, term by term. Using integration 
by parts in time we note that 




{A ^utt){h-Vu) 


(d^(A ^ut),d„u)j. 


T 

0 



{h-v){d„{A ^ut){d^ut). 


From the elliptic regularity theorem [531 Ch. 3]: for any h G T^(U), we note ||A ^h \\2 < Ull/iHo.n, 
and thus we have: 


\mA-^h)\\o,r< C\\A-^h\\y2+e,n< C\\A-^h\\2,a< C\\h\\o,a. 


Hence, compactness of the Sobolev embeddings (and Lions’ Lemma [501 P-108]), the Holder-Young 
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inequality yields that 


/ {d,{A-^Uu){h-yu)<eEiT)+C{e)\\u{t)\\o,n+CE{0)+e i?(r)+C(e) / . (4.17) 

Jo Jr Jo Jo 

Again, from standard elliptic theory, di,Di^ £ ££ (L^(r), i7“^(r)). Thus we have 

J j (Am|p) (h • Vu) <J ||h • Vu|| 


ffi(r) 


< 


c(h)(||zi(a,u*)ll L2(0,T;L2(r))) ^ 


d'^i 


dv"^ 


L\T) 


d\ ' 


drdv 


L^T) 


1 


1/2 


and similarly, 


[ [ d^{A ^p(x))(/i • Vu) < e / E{t) + C{p,e)-T. 

Jo Jr Jo 


(4.18) 


(4.19) 


Taking into account (I4.17|) - (I4.19D in (14.161) . considering u = 0 on T, employing Assumption [U using 
dm, and applying the Holder-Young we thus arrive at the next preliminary estimate which can be 
implemented in (14.11^ : 


-2 


f ( d^{Au){h -Vu) - [ ( 7 -||Vu||^) f {h ■ iy){d^u)'^ 
Jo Jr Jo Jr 


(4.20) 


/o Jr 

< 


C{h) {\\Did,ut)\\ L2(0,T;L2(r))) 


d'^u 


dv'^ 


L 2 (r) 


d\ ' 


drdv 


1/2 


L 2 (r)_ 


+ C(e)y^ \\d,ut\\l,^^^+e^E{T) + E{t)'^ + CE{Q) + C{p,e,-i,h)-T 

+ 7 / [ {h ■ 

Jo Jr 


(4.21) 


For the terms in the last line of (14.111) involving the higher order trace term d^{h ■ Vit), we note the 
boundary condition Am = —D{d^ut) on T , and we have 


/o Jr 


D{d^ut){d^u) 


lo Jr 


2{Au)d„{h ■ Vu) < 


\\D{dvUt)\\L^(o,T-,r)\\duu\\L2(^Q j,.j.^ + \\D{duUt)\\L^{o,T-L^{r))\\du{h ■ VM)||i2(o_y.^2(pp. (4.22) 
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We now combine ()4.21|1 and (j4.22|) in (I4.11|l . and absorb terms to obtain: 



E{t) < C{h) 



d'^u 


2 

+ 

L 2 (r) 


d'^u 

drdv 


2 

L 2 (r) 


1/2 


+ C(e) r \\d,ut I|i 2 (r) + e^(T) + CE{0) + C{p, e, 7 , h)T 
Jo 

+ C\\D{d,jUt)\\L^(o,T-,r)\\duu\\L2(o,T-,r) + C\\D{dvUt)\\L'^(o,T-,L^(r))\\dv{h ■ VM)||z,2(o_T;z,2(r)) 

(4.23) 


Now, considering (14.111) over the interval (a, T — a) instead of (0, T) (hence, performing the calculations 
above on (a,r — a)), we can use the decreasing nature (modulo a constant) of the energy functional 
( (|1.10l) and (14.31) ) . and then extend the non-critical integrals back on (0,T). We obtain: 



E{t) < C{h) ^||.D(9[/Ut)||p2(Q 7 ’.p 2 (p))^ 

rT 



d'^u 


2 

L 2 (p) 


d'^u 

drdv 


L 2 (p )_ 


+ C(e) / \\d,jUt\\'^ 2 fT^\ + €E(T — a) + CE{a) + C{p,e,j,h) ■ [T — 2a] 

Jo 

+ C\\D(d,yUt)\\L^(o^T-,r)\\dvU\\L^(a,T-a-r) 

+ C\\D(d,yUt)\\L^(o^T-,L^(r))\\du{h ■ Vm)||l 2 (q,_p_„.p 2 (p)) 


1/2 


(4.24) 


The following sharp regularity result for the boundary traces of solutions to the Euler-Bernoulli 
equation (linear, with given RHS) will be critically used: 

Theorem 4.3. fflSl p.281,Theorem 2.4 ]) Let Q<a<T,Q<5< 1/2, and sq < 1/2 be arbitrary. 
Then the generalized solutions of cm with (HD) boundary conditions enjoy the bound 


/ ^-oc n 1/2 , 

J \drru\'^+ \d,„,u\'^+ \dr,.u\‘^drdt^ < C'(T,a,( 5 )|||/(u)||p 2 ( 0 _p.^~«o(n)) 

+ ll^llL2(0,T;if2-5(Q)) + ||Zl(9j^Mt)||p2(o,T;L2(r)) + ll'9j^^^t||L2(o,T;L2(r))| 


(4.25) 


Hereafter, the explicit dependence of the above constants on a, S, 7 , h will be suppressed. For the 
first term of RHS of above inequality, we make use of (14.251) and recall the nonlinear energy E] from 
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Young’s inequality we obtain: 


||-D(9,.Ut)IU2(0,T;L2(r))| J 


d'^i 


dv^ 


+ 


L=(r) 


d'^i 


drdy 


1/2 


L=(r) 


<C + C(T)|p(a.uO||i.(o,T;L^(r)) / ||(7-||Vu|nAu||:^^dT 


■-Jo 


C{T) ||£'(i9i,Mt)||^2(o,T;L2(r)) + ll^i"'“t|lL2(0,T;L2(r)) + II“IIl2(0,T;H2-«(o)) 


p'^ 2 

<C + C(T)|p(5.uO||i.(o,T;L^(r))[y^ (^(r)) dr 


C{T) ||£^(i9jy'Ut)llL2(0,T;L2(r)) + I + C'l |w| |p2(Q (q)) 


rT-c 


<C + e E{t) + CE{a) + C{T,e) + C{T)\\D{d,^ut)\\'j^ 2 (^Qj-,L^{r)) 

J cx. 

+ C'(T) + ll^i^'^i|lL2(;o,T;L2(r)) 


■-^0 


(^(l 


r) 1 dr 


(4.26) 


We proceed with estimating the last line of the RHS of (I4.24p . Using Young’s inequality, the trace 
theorem, and (HU we get 


pT — Ol pi 

||di>(i9j,Ut)||L2(o^T;r)||9i/u||i2(„^p_Q,.r) < e / E{T)+CE{a)+ / \\D{dvUt)\\\' 2 ^Y)+C{e)-T. 

J OL J 0 

Considering 

dv{h ■ Vm) = dvU + {h ■ i')d„yU + {h ■ T)d^rU, (4.27) 

and applying similar steps as before, we have: 

\\D{di^ut)\\L^(o,T-,L^{r))\\dv{h ■ Vu)||p2(„_p_Q,.p2(r)) < 

e f E{t) + CE{a) + C{T) + C{T)\\D{d^ut)\\\2(^Qrp.^2(^Y')) 

J CX. 

+ C{T) ||7?(9jyUt)|||2(o,T;L2(r)) + ll'9i'Wt|||2(o_T;L2(r)) 

If we take into account the last two inequalities we have: 

||£)(9j/Ut)||L2(o,T;r)||c^!2u||i2(„_p_Q,.r) + \\D{dvUt)\\L‘^(o^T-L^(T))\\dv{h ■ Vu)||l2(q, p_„.p2(r)) 

<e r “.g(r) + CU(a) + C(T) + C(T)|p(5.u0lli2(o.T;L2(r))[ {E{T)fdr 


'-Jo 


{my 


dr 


'-JO 


+ C{T) \\D{di,ut)\\‘j^2(^Q x-L'^(r)) + ll'^‘2Ut|li2(Q p.p2(pp 
Finally, we recall that from the energy identity and (ll.lOp . we have for s < t: 

E{t) < C'iU(s) + C 2 
[t — s]U(f) < Cl f U(r)dr + C 2 [t — s] 

J S 


(4.28) 

(4.29) 

(4.30) 
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Thus, integrating the energy equality (I1.10|) on [a,T — a] and employing (j4.3|) . as well as using ()4.26|1 
(j4.28|l in (j4.24L we finally have: 


(T - 2a)E{T) < CE{Q) + / {E{T)fdT 


'-Jo 


+ C'(T) ||T'(9,yMt)|||2(o,T;L2(r)) + l|l9i'Ut|||2(o,T;L2(r)) +C{T), 


which is the desired observability estimate, written in terms of the nonlinear, positive energy E. □ 

Remark 4.4. Going back to the equation to obtain trace estimates (and making use of the sharp trace 
result (14.251) above) is the key insight into the present problem. This step connects our observability 
type estimate to the long-time behavior analysis of m by utilizing the fact that the nonlinear energy 
appears under the time integration against the damping mechanism. 

Remark 4.5. We emphasize that at multiple times in the argument above we rely critically on the 
structure of the Berger nonlinearity —particularly in the ability to drop the term 

- f ||Vu|p [ {h ■ 

Jo Jr 

via the geometric condition and in the estimate of higher order trace terms (See Remark |4.41) . Since 
the von Karman nonlinearity does not possess a structure which accommodates the trace estimates 
needed above, we cannot expect a straightforward approach to yield similar results there. Thus, for 
the Berger nonlinearity, in addition to the existence of a compact global attractor, we also have an 
explicit estimate on the absorbing set (not available when the attractor is obtained—for gradient 
systems—indirectly through the asymptotic smoothness property). 


Completion of the Proof of Theorem \4-l\ 

Now, let us denote 


Dq = f [ D{duUt){dyUt). 
Jo Jr 

Since D{0) = 0, we have for s > 2 that 


(4.31) 


D'(Ts)dT > sinf [ D'(Ts)dT 

^>■2 Jo 


> sinf [ D'(Ts)dT > -inf inf D'(ts) = -infZ?'(£) > —m. 

s>2 ji /2 ^ ^ 2 s >2 1 / 2 < t <1 ^ 2C>1 ^ 2 

It can be shown in an analogous way that for s < —2 we have D{s) < |to. Additionally, with a given 
e > 0 we have for s G (e, 2): 



D{s) 


f D'{T)dT > f D'{T)dT > me 
Jo Jo 


s 

2 
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with a same inequality for s € (—2, —e). Hence the above relations give that: 


D(s)s > —m for |s| > 2 


and D{s)s > rue 


for |s| > e . 


Similarly, one can show for s > 0 that 


D{s) < s ■ max 


sup D\^),M 
0<|?|<1 


Hence we obtain 


{D{s)) < D(s)smax 


sup 

0<|«|<1 


M 


Now, considering the above relations we obtain that 




T;L2(r )) 




T;L2(r )) 


< C'(T, r) + max ■ 


, sup 
0<|«|<1 


M)dI (4.32) 


If we take into account (|4.32l) in our observability estimate (14.51) we have 


{T -2a)E{T)dT <CE{Q) + C{T)Dl+ C[T)Dl / {E{T)f + K{T) 


where K{T) does not depend on ii^(O). Since 


E{Q)<Dl + CE{T) + C{p,^), 


we have 


(T - 2a - C)E{T) < C{T)D^ 



+ K{T). 


Now, if we rewrite the last inequality in terms of the full nonlinear enery iS‘{T), use the relation 
(1431) between S’ and E, and employ the notation Sm{T) = S{T) + M, where M is the constant coming 
from (14.41) . we obtain: 


{T-2a-C)S m{T)<C{T)dI 



+ K{T) 


Using the fact that Sm{T) < ^m(O), and Dq = Sm{0) — Sm{T), we have 


{T-2a- C)Sm{T) < C{T) [Sm{0) - Af(r)] [l + TS^{0)] + K{T), 
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from which we obtain 


(T - 2a - C) + C{T) [1 + T^^(O)] <gM{T) < C{T) [l + r^^(O)] <fM(0) + K{T). 


Then, assuming T large enough and recalling (14.3|) . we have 


T + C{E{Q),T) 


T + C{E{Q),T) 


Denoting r]{E{0),T) = < 1, we rewrite the last inequality as 


^m{T)<v{E{0),T)^m{0) 


K{T) 


T + C{E{0),T) 


Now, we can reiterate the same estimate on each subinterval {mT, (m+l)T) via the semigroup property. 
We note that the constants C{E{0),T) and K{T) will be same at each step. Then we obtain 


a™ , ^ C{E{ 0 ),T) 

&M{{m + Iji ) <-- &M{rnl ) 


KIT) 


T + C[E{Q),T) T + C[E{Q),T) 

< v{E{0), T)-A^(O) + f2 viE{0),T)^ 

1 


i=0 


T + C{E{0),T) 
K{T) 


<v{E{Q),T)”^Sm{Q)- 

- r]{E{0),T)^ T + C{E{0),T) 

= 77(i?(0),T)™<r(0) + i?(T) 

The monotonicity and continuity of the energy function, and the fact that r]{E{0),T) < 1, gives 

<^(T) + 1, for all t>to{E{0)), 
which together with (j4.3|) completes the proof of Theorem 14.11 


4.2.3 Asymptotic Smoothness—Completion of the Proof of Theorem 11.31 

In this section, as the second step of the proof of Theorem 11.31 we show the asymptotic smoothness 
property of the dynamical system generated by HU with (HD) boundary conditions in the following 
theorem: 

Theorem 4.4. Let Assumptions[Ji and\^ be in force. The dynamical system generated by 

generalized solutions to HU taken with (HD) boundary conditions is asymptotically smooth. 

We note that after Theorem 14.41 is proved, our main theorem on the existence of a global attractor 
for generalized solutions to (HU taken with (HD) conditions lTheorem ll.3|l will be proved, via Theorem 

O 

Proof of Theorem \4.4\ The proof rests on the application of Theorem 16.31 Hence we are interested 
in the difference of two solutions z = u — w, where U{t) = {u{t),ut{t)) = SH{t)yi and W{t) = 
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{w(t),wt{t)) = SH(t)y 2 solve (jl.lll corresponding to initial conditions yi and y 2 (respectively), taken 
from an invariant, bounded set. Then, z will solve the following problem: 


Ztt + A^z + J^(z) = 0 in(0,T)xff; z = 0, IS.Z =—[D{dvUt) — D{di,wt)\ on F, (4.33) 
z(0)=uo-w;o, Zt{Q)=ui-wi, 

where T{z) = /b(u) —/b(w). We denote the energy for the difference: E^{t) = i{||Az|p + ||zt(t)|p}. 
The energy identity (11.101) . the relation (14.31) between the energies, and the bound on lower frequencies 
gH, give that there exists an i?* such that the set 


= {(mo,ui) G : S{uo,ui)<R} 

is a non-empty bounded set in for all R> R^,- Moreover, any bounded set ^ C is contained in 
for some R, and the set Wr is invariant with respect to Snit). Then, we consider the restriction 
of the dynamical system {Sh{-),^) to {Sh{-),'^r) in showing the asymptotic smoothness property, 
and thus we consider the solutions u, w satisfying 

\\u{t )\\2 + l|ut(i)||o + \\w{t )\\2 + ||wt(t)||o < C{R), t > 0. 


The main ingredient will be the following estimate on the difference of the solutions: 

Lemma 4.5. Let T > 0 and /3 G C^(]R) be a given function satisfying (i) supp{/3) C [a,T — a] (with 
a < Tf2); (ii) 0 < /3 < 1 and /3 = 1 on [a,T — a]. Then, any solution z to equation (14.331) satisfies 


ES)P{t)dt < Cl / E,{t)\P'{t)\dt + C2{T) 


D{d„ut) - D{d,,wt) d,,Zt + \d,,Zt\'^dT 


/o Jr ■- 


-|-/.o.t.(z). 


(4.34) 


where 


l.o.tfz) = C{T,R) sup ||z(t)|| 2 _,, + 

L [O.T] 


l^t('r)ll%dr , 0<77<1/2. 


Proof. The proof of this lemma follows m Section 10.5] exactly, with respect to the linear portion 
of the dynamics and the nonlinear damping: the application of the flux multiplier h ■ V(/lz), where 
h = X — xq, Xo G the use of sharp trace results proved in (as in the previous section), the 
local Lipschitz property of /, and the employment of Assumption [TJ This yields the above estimate 
(14341) . □ 


Now, by the energy relation, 

rT 

/r 

we have: 


E.{T)+ f f [D{d,.ut) - Did^^wt)] {d,.zt) + [ iEiz),Zt) = E,{t), T>t, (4.35) 
J t Jr J t 
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(i) 


ES){i-m + wmdt< 


where 


(ii) 


[l-Pit) + \P'it)\]dt E,{T)+ / [Did,ut)-D{d,wt)]id,zt) +Mz), 

^ Jo Jt 

E 4 z)= [ ( 1 -/ 3 (i) - |/ 3 '(i)|) [ {E{z),zt)dTdt 
Jo Jt 

TEz{T)< [ Ez{t)dt-\- 

Jo 


where 


E^^{z) = - [ [ {E(z),Zt)dTdt. 

Jo Jt 


Combining (i) and (ii), we obtain that there exist Tq > 0, and constants Ci{T) and C2{R,T), such 
that 


TE,{T)+ f E,{t)dt < Ci{T) [ [ [D{d,ut)-D{d,wt)]{d,zt)dr dt+ [ [ \d,ztf dV dt 

Jo Jo Jr Jo Jr 

+ E^{z)+J^^^{z) + C 2 {R,T)l.o.t.{z), V T>To. (4.36) 

Using the assumption on the structure of the damping (Assumption [TJ , we have 

[ [ \dyZt\^ <e^C{e) [ ( [D{d^ut) - D{d^wt)\{dyZt), V e > 0. 

Jo Jt Jo Jt 

On the other hand, it follows immediately from the energy relation (14.351) that 

[ [ [D{d,Ut)-D{d,wt)]{d,Zt)<EM-E,{T)+ [ {E{z),zt) 

Jo Jr Jo 

Now, taking into account the last two inequalities in (14.3611 . we obtain 

TE,{T)+ [ E,{t)dt<e + Ciie,T)[E,{0)-E,{T)] + C 2 {e,T) [ (E(z),Zt) 

Jo Jo 

+ J^4z)+E^4z)+C3{e,T,R)Lo.t.iz) (4.37) 


for any e > 0 . 

In order to apply Theorem 16.31 we need to simplify the nonlinear terms above. For this, we 
demonstrate a decomposition of Berger nonlinearity for the difference of two solutions. 

Lemma 4.6. Let z = u — w, and let /b(u) = (7 — ||Vu|p)Au, and E(z) = /b('«) — /b(w). Also 
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assume u,w G C{[t,T]] (H^ fl i^Q )(fJ)) fl C^{[t, T]-,L^{fl)). Then 

r{T{z),Zthdr=[iT{z),z)n + ^WVzf - l\\Vu\f\\yz\\^ - (||V«||2 - ||Vzi;||2)(Azz;,z)n 


1 ^ 
- t 


- ( (||Vu||^ - ||Vw||^)(Aw,zt)ndr + f {Au,ut)n\\Vz\f dr 

Jt Jt 


Proof. Letting z = u — w, we note 


[ {T{z),Zt)dT = - f ([||Vm||^Am - ||Vw;||^AM;],2:t)dr+ 7 / {Az,Zt)dT 

Jt Jt Jt 

= [ ([||Vm|PAu]'— [||Vu>|pAii;]', z)dr + {T[z), z)-\-'^\\\/z\ 


where we have integrated by parts in time, integrated by parts in space on the term with 7 , and 
recognized the total derivative term (Vz, Vzt). We now focus on the nonlinear term and observe 

([||Vw|PAm]' — [||Vr(;|pAu>]',z)dr ((||Vu|p)'Au — (||Vw|p)'Ar(;,z)dr 

If we add and subtract the terms which are mixed in u and w we have 
f ([||Vu|pAu]' — IWVwW^Aw]', z)dT = 


rT n 


(||Vu|n'(Az,^)+[(||Vu||7 - (livin'] (Arc, z) 


dr 


rT , 


WVuW^iAzt.z) + [||Vm|P - ||Vw|p](Awt,z) 


dr. 


Now, recognizing a total derivative and adjusting we arrive at 

pT 


([llVulpAu]' — [||Vw|pAw]', z)d7 


[^[l|Vu|n|V^||V^I|V«ll"(l|Vz||2)' 


dr 


rT ^ 


ft 

rl 


d 

-dr 


[(||Vu||2 - ||Vu;||2)(Au;,z)] + (||Vu||2 - ||Vu;||2)(Au;,z*) 

T 

t 


dr 


= - _2l|V^II1|V^lr + (llVull^ - ||Vu;||^)(Au;,z) 


f (||Vu||^ - ||Vw||^)(Aw,zt)dr + / (Au,Mt)||Vz||^dr 

Jt Jt 


which is the desired decomposition. 


□ 
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At this point, noting that 


||m(<)I |2 + l|ut(i)||o + ||w(0ll2 + lkt(i)||o < C{R), t > 0, 


using triangle inequality 


||Vz.||2-||Vn;|n = ||Vu||-||VHI (||Vu|| + ||Vn;||) 
< ||Vm-Vw||(||Vu|| + ||Vw||) 
<C(i?)||z||i, 


and taking into account the last two inequalities in Lemma 14.61 we obtain the key inequality: 


{R{z),Zt}dT 


< e 


[ Ez{T)dT+ C{e,R,T) sup\\z\\l_ , r] > 0. 
Jt [t.T] 


(4.38) 


From this, the analogous bounds on and follow. 

Completion of the Proof of Theorem \4.4\ ^ Asymptotic Smoothness 

We are now in a position to finish the proof of Theorem 14.41 Firstly, implementing (|4.38ll in (I4.37|) 
and rescaling e, we see that 


E,{T)<e + Ci{e,T)[E,{0)-E,{T)] + l.o.t.{z). (4.39) 

In order to invoke Theorem l6.3l we need only to construct a functional d'e^ satisfying the compensated 
compactness condition [HEI], i-e., 

liminf liminf '^e,R,T{yn,yk) = 0 

k—fcx) n—^oc 


for every sequence {yn} from C Jif. For this, we use a standard stabilization type argument via 
the semigroup property of the dynamics on the interval ((m — 1)T, mT). This yields 

EzimT) < ^ ^ ^ + j,E,{{m - 1)T) + l.o.t.^iz), 

Cie,R,T) 


where 


le = 


l + C{e,R,T) 


< 1 


and 


l.o.t."^{z) = C{€,R,T) sup \\z(t)\\1_^ + 

*- [(m —l)T,mT] 


pmT 


\\zt{T)\\t^dT 


After iteration, we observe that 


m —1 


E^imT) < e + 7r^^.(0) + ^ (z). 

j=o 
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For any eg, we may select e > 0, then m large enough so that, with H and T fixed, we have 


m —1 


Ez{mT) < eo + ^ yf/.oi.™ ^{z). 
3=0 


Now, for 2/1 = [/(O) and 2/2 = hF(0), setting 

m—1 

7^/.O.t (z) 

3=0 

the functional 4'eo,^,mT—being compact with respect to the energy space (via the Sobolev embeddings)— 

satisfies the requisite condition of the functional 'h in Theorem 16.31 This concludes the proof of 
asymptotic smoothness of the dynamical system (5'ij(•), #){) for any i? > 0, according to Theorem 

[Ql □ 

As a result, by Theorem 14.11 and Theorem 14.41 we obtain the existence of global attractor which 
finishes the proof of Theorem 1 1.31 
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6 Appendix 

6.1 Long-time Behavior of Dynamical Systems 

In the context of this paper we will use a few keys theorems (which we now formally state) to prove 
the existence of the attractor and determine its properties, as well as provide some context for other 
results mentioned in the discussions above. For general dynamical systems references, see [inTiisniiii] 

(and refences therein). For proofs pertinent to what is presented here, and more references, see |14| . 

Let {H, St) be a dynamical system on a complete metric space Ji with ^ = {x G H : Stx = 

X for all t > 0} the set of its stationary points. (Ji, St) is said to be dissipative iff it possesses a 
bounded absorbing set B. This is to say that for any bounded set B, there is a time ts so that 
StsiB) C B. We say that a dynamical system is asymptotically compact if there exists a compact set 
K which is uniformly attracting: for any bounded set D (ZTi we have that lim d-u{StD\K} = 0 in 
the sense of the Hausdorff seinidistance. (H, St) is said to be asymptotically smooth if for any bounded, 
forward invariant {t > 0) set D there exists a compact set K <Z D which is uniformly attracting (as 
in the previous definition). Global attractor A is a closed, bounded set in 'H which is invariant (i.e. 

StA = A for all t > 0) and uniformly attracting. 

The following if and only if characterization of global attractors is well-known mm 
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Theorem 6.1. Let (H,5't) he a dissipative dynamical system in a complete metric space TL. Then 
{TL,St) possesses a compact global attractor A if and only if{'H,St) is asymptotically smooth. 

An asymptotically smooth dynamical system for which there is a Lyapunov function $( 1 ) that is 
bounded from above on any bounded set can be thought of as one which possesses local attractors. To 
see this stated precisely see m page 33. Such a result provides an existence of local attractors, i.e., 
and attractor for any bounded set of initial data. However, these sets need not be uniformly bounded 
with respect to the size of the set of initial data. The latter can be guaranteed by the existence of 
a uniform absorbing set. However, establishing this existence of an absorbing set may be technically 
demanding. In some instances, there is a way of circumventing this difficulty which takes advantage 
of the “good” structure of a Lyapunov function. 

A strict Lyapunov function for {fH,St) is a functional <i> on 'H such that (i) the map t —>■ ^{Stx) 
is nonincreasing for all a; € "H, and (ii) ^{Stx) = $(a:) for all t > 0 and x € TL implies that a; is a 
stationary point of {T-i,St). If the dynamical system has a strict Lyapunov function defined on the 
entire phase space, then we say that (fH, St) is gradient. 

We can address attractors for gradient systems and characterize the attracting set. The following 
result follows from Theorem 2.28 and Corollary 2.29 in m- 

Theorem 6.2. Suppose that {'H,St) is a gradient, asymptotically smooth dynamical system. Suppose 
its Lyapunov function $(x) is bounded from above on any bounded subset ofH and the set = {x G 
TL : <i>(x) < R} is hounded for every R. If the set of stationary points for (T-L,St) is bounded, then 
{TLjSt) possesses a compact global attractor A which coincides with the unstable manifold, i.e. 

A = = {x GH : 3 [/(t) S "H, V t G M such that U{0) = x and lim d-ulJJ{t)\^) = 0}. 

t—>■ —00 

Secondly, we state a useful criterion (first appearing |21) and stated in the present version in 
M) which reduces asymptotic smoothness to finding a suitable functional on the state space with a 
compensated compactness condition: 

Theorem 6.3. f\13^ -Proposition 2.10) Let ['H.,St) be a dynamical system, TL Banach with norm || • ||. 
Assume that for any bounded positively invariant set B G TL and for all e > 0 there exists a T = 
such that 

IlS'rXi - STX2\\n < e + e,B,T{xi,X2), Xi € B 

with a functional defined on B x B depending on e, T, and B such that 

liminf liminf = 0 

m n 

for every sequence {xn} C B. Then (TL,St) is an asymptotically smooth dynamical system. 

References 

[ 1 ] G. Avalos, P.G. Geredeli, and J.T. Webster, Finite dimensional smooth attractors for the Berger 
plate with dissipation acting on a portion of the boundary, preprint. 


34 


[2] G. Avalos and 1. Lasiecka, Exponential Stability of a Thermoelastic System Without Mechanical 
Dissipation, Rend. Istit. Mat. Univ. Trieste, 28 (1997), pp. 1-28. 

[3] G. Avalos and I. Lasiecka, Boundary Gontrollability of Thermoelastic Plates via the Free Boundary 
Gonditions, SIAM J. Control. Optim., 38 (2000), pp. 337-383. 

[4] A. Babin and M. Vishik, Attractors of Evolution Equations, North-Holland, Amsterdam, 1992. 

[5] B. Banerjee and S. Datta, A new approach to an analysis of large deflections of thin elastic plates, 
J. Non-Lin. Mech., 16 (1981), 1, pp. 47-52. 

[6] H.M. Berger, A new approach to the analysis of large deflections of plates, J. Appl. Mech., 22 
(1955), pp. 465-472. 

[7] L. Bociu and D. Toundykov, Attractors for non-dissipative irrotational von Karman plates with 
boundary damping, .IDE, 253 (2012), 12, pp. 3568-3609. 

[8] V.V. Bolotin, Nonconservative Problems of Elastic Stability, Pergamon Press, Oxford, 1963. 

[9] S.G. Brenner and R. Scott, The mathematical theory of hnite element methods, v. 15, Springer 
Science & Business Media, 2008. 

[10] I. Chueshov, Introduction to the Theory of Infinite Dimensional Dissipative Systems, Acta, 
Kharkov, 1999, in Russian; English translation: Acta, Kharkov, 2002; 

see also http://www.emis.de/monographs/Chueshov/ 

[11] I. Ghueshov, Introduction to the Theory of Infinite-Dimensional Dissipative Systems, Acta, 
Kharkov, 1999 (in Russian); English translation: Acta, Kharkov, 2002; 

http://www.emis.de/monographs/Chueshov/ 

[12] I. Chueshov, M. Eller, and I. Lasiecka, Finite dimensionality of the attractor for a semilinear wave 
equation with nonlinear boundary dissipation. Comm. PDE 29 (2004), pp. 1847-1976. 

[13] 1. Chueshov and 1. Lasiecka, Long-time behavior of second-order evolutions with nonlinear damp¬ 
ing, Memoires of AMS, v. 195, 2008. 

[14] 1. Chueshov and 1. Lasiecka, Von Karman Evolution Equations, Springer-Verlag, 2010 

[15] P. Ciarlet and P. Rabier, Les Equations de Von Karman, Springer, 1980. 

[16] A. Eden and A.J. Milani, Exponential attractors for extensible beam equations. Nonlinearity 6 
(1993), pp. 457-479. 

[17] P.G. Geredeli, 1. Lasiecka, J.T. Webster, Smooth attractors of finite dimension for von Karman 
evolutions with nonlinear frictional damping localized in a boundary layer, J. Diff. Eqs. 254 (2013), 
pp. 1193-1229. 

[18] P.G. Geredeli, J.T. Webster, Decay rates to eqilibrium for nonlinear plate equations with geo¬ 
metrically constrained, degenerate dissipation, Appl. Math, and Optim. 68 (2013), pp. 361-390. 
Corrigendum, Appl. Math, and Optim. 70 (2014), pp. 565-566. 


35 



[19] G. Ji and I. Lasiecka, Nonlinear Boundary Feedback Stabilization for a Semilinear Kirchhoff Plate 
with Dissipation Acting Only via Moments-Linhting Behavior, JMAA 229 (1999), pp. 452-479. 

[20] S. Kesavan, Topics in Functional Analysis and Applications, John Wiley and Sons, India, 1989. 

[21] A.Kh. Khanmamedov, Global attractors for von Karman equations with non-linear dissipation, 
J. Math. Anal. Appl 318 (2006), pp. 92-101. 

[22] H. Koch and I. Lasiecka, Hadamard well-posedness of weak solutions in nonlinear dynamic 
elasticity—full von Karman systems. Evolution Equations, Semigroup and Functional Analysis, 
50 (2002), pp. 197-212. 

[23] J. Lagnese, Boundary Stabilization of Thin Plates, SIAM, 1989. 

[24] I. Lasiecka and R. Triggiani, Control Theory for Partial Differential Equations, Cambridge Uni¬ 
versity Press, Cambridge, 2000. 

[25] 1. Lasiecka and R. Triggiani, Sharp Trace Estimates of Solutions to Kirchhoff and Euler-Bernoulli 
Equations, Appl. Math Optim 28 (1993), pp. 277-306. 

[26] T.F. Ma, Boundary stabilization for a non-linear beam on elastic bearings. Math. Meth. Appl. 
Sci. 24 (2001), pp. 583-594. 

[27] A. Miranville, S. Zelik, Attractors for dissipative partial differential equations in bounded and un¬ 
bounded domains. In: M.C. Dafermos and M. Pokorny, (eds) Handbook of Differential Equations: 
Evolutionary Equations, v.4, Elsevier, Amsterdam, 2008. 

[28] J.L. Nowinski and H. Ohnabe, On certain inconsistencies in Berger equations for large deflections 
of plastic plates, Int. J. Mech. Sci 14 (1972), pp. 165-170. 

[29] J.-P. Puel and M. Tucsnak, Boundary stabilization for the von Karman equations, SIAM J. 
Control and Optim. 33 (1995), 1, pp. 255-273. 

[30] G. Raugel, Global attractors in partial differential equations. In: Fiedler, B. (ed.) Handbook of 
Dynamical Systems, v. 2, Elsevier Sciences, Amsterdam, 2002. 

[31] R.E. Showalter, Monotone Operators in Banach Spaces and Nonlinear Partial Differential Equa¬ 
tions, V. 49, AMS, 1997. 

[32] C.P. Vendhan, A study of Berger equations applied to nonlinear vibrations of elastic plates, Int. 
J. Mech. Sci 17 (1975), pp. 461-468. 


36 


